Optimizing Anharmonicity in Nanoscale Weak Link Josephson Junction Oscillators 
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Josephson tunnel junctions are widely used as nonlinear elements in superconducting circuits 
such as low noise amplifiers and quantum bits. However, microscopic defects in the oxide tunnel 
barrier can produce low and high frequency noise which can potentially limit the coherence times 
and quality factors of resonant circuits. Weak link Josephson junctions are an attractive alternative 
provided that sufficient nonlinearity can be engineered. We compute the current phase relation 
for superconducting weak links, with dimensions comparable to the zero temperature coherence 
length, connected to two and three dimensional superconducting electrodes. Our results indicate 
that 50-100 nm long aluminum nanobridges connected with three dimensional banks can be used to 
construct nonlinear oscillators for bifurcation amplification. We also show that under static current 
bias, these oscillators have a sufficiently anharmonic energy level structure to form a qubit. Such 
weak link junctions thus present a practical new route for realizing sensitive quantum circuits. 



Nonlinear oscillators with low loss are a basic building 
block of superconducting quantum information process- 
ing circuits. They are used both as "artificial atoms" 
to realize quantum bits and as the gain element 
in low noise amplifiers Q for quantum state measure- 
ment. The Josephson tunnel junction, formed by sep- 
arating two superconducting electrodes with a barrier 
such as a thin oxide layer [3j, is the nonlinear circuit 
element commonly used to construct these anharmonic 
oscillators. The current I(t) and voltage V(t) of the 
junction can be expressed in terms of S(t), the gauge- 
invariant phase difference, as I(t) = IosmS(t) and 
V(t) — (po d8/dt, where the parameter Iq is the junc- 
tion critical current and tpo = h/2e is the reduced flux 
quantum. These equations parameterize a nonlinear in- 
ductance Lj — tpo/ {Iq cos(<5)}. This Josephson induc- 
tance is shunted in parallel by a capacitance Cs resulting 
from the geometric capacitance of the junction and possi- 
bly an additional external capacitor, thereby forming an 
electrical oscillator with an amplitude dependent plasma 
oscillation frequency up = 1/ \J ' L,j{I)Cs- In principle, 
Josephson oscillators should exhibit a very large qual- 
ity factor (Q) since at temperatures T well below the 
critical temperature T c , dissipative thermal quasiparti- 
cle production is exponentially suppressed. In practice, 
microscopic defects in the tunnel barrier 4| and in the 
shunt capacitor dielectric [B[ can reduce Q. 

An alternative structure to the oxide barrier tunnel 
junction is a weak link Josephson junction (see [6| for a 
review) where a small metallic constriction with dimen- 
sions of order the superconducting coherence length £(T) 
bridges two large superconducting electrodes. This type 
of junction also behaves as a nondissipative inductance 
but the strength of the nonlinearity is a sensitive function 
of the weak link and electrode geometry. In contrast to a 
tunnel junction which invariably has a sinusoidal current 
phase relation (CPR) with / oc sin(<5), a weak link junc- 



tion at T << T c can exhibit a wide range of CPRs rang- 
ing from I oc S — characteristic of a linear inductor — to 
a distorted sinusoid I oc cos (5/2)tanh" 1 [sin( l 5/2)] for an 
idealized, short constriction [7j. In this letter, we present 
calculated CPRs, obtained via a numerical solution of the 
Usadel equations [8[, for thin, diffusive superconducting 
bridges contacted with two and three dimensional banks 
for varying bridge length and width. We also compute 
the dynamics of a weak link Josephson oscillator under ac 
excitation, and its bound state energies. Our results in- 
dicate that with a 50-100nm long, 45nm wide aluminum 
bridge with three dimensional banks, the oscillator is suf- 
ficiently nonlinear to observe a stable dynamical bifurca- 
tion with a fixed frequency microwave drive, and to ob- 
tain three quantized levels with ~ 3% anharmonicity un- 
der dc bias. These results are obtained without putting 
impractical constraints on the shunting capacitance, os- 
cillator quality factor, or current bias stability, and thus 
suggest a new route for realizing superconducting qubits 
and amplifiers without tunnel junctions. Furthermore, 
this junction geometry is suitable for coupling Joseph- 
son devices to other solid state quantum systems such as 
nanomagnets and quantum dots. 

We consider two canonical weak link geometries: a 
"Dayem" bridge where the banks are the same thickness 
as the weak link and a "variable thickness bridge" where 
the banks are significantly thicker than the weak link. 
Considerable literature exists on both of these struc- 
tures. As T — > T c , the temperature dependent coher- 
ence length £(T) diverges and mean field methods, such 
as Ginzburg-Landau theory Q, can be used to compute 
the magnitude and phase of the superconducting order 
parameter along the bridge and in the banks [1C| . For 
qubits and low noise amplifiers, we are interested in op- 
crating the Josephson junction at T <C T c to minimize 
quasiparticle loss. In this regime, linear response meth- 
ods based on Gorkov functions have been applied to short 
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FIG. 1: Computed current phase relations for an Al 
nanobridge with width W — 45 nm and three different lengths 
L — 15, 45, and 240 nm. Data are shown for (a) 2D and (b) 
3D banks at T = 0.15 K. The 3D structures exhibit nonlinear 
CPRs which progressively become linear as L exceeds £. Cor- 
responding lengths for 2D structures exhibit a nearly linear, 
hysteretic CPR (the other branch is suppressed for clarity). 
In both geometries, the maximum current systematically de- 
creases with increasing bridge length. 



weak links with length L and width W <C £ (T) , and with 
rigid boundary conditions imposed at the bridge ends — 
Kulik-Omelyanchuk (K-O) theory Q. However, 10 — 30 
nm thick Al bridges with a mean free path I ~ 1 nm flT| 
have £(0) » 30 nm, making the short limit difficult to 
achieve with conventional e-beam lithography. We thus 
focus on weak links with lateral dimensions comparable 
to £(T = 0) and large compared to I. In the absence of 
closed form expressions for the CPR in this parameter 
regime, we apply the Nambu-Gorkov formalism [l2|, HH 
in the diffusive limit and numerically solve the Usadel 
equations [8j for two and three dimensional weak link 
junctions. 

We first compute the complex pairing potential A in 
the bridge and the banks. The modulus and argument of 
A yield the superconducting gap and phase, respectively. 
We use the $ parametrization of the Usadel equations 
1141 by making the substitutions F = —, — 2 and G 2 = 

1 — |F| 2 where G and F are the normal and anomalous 
Green's functions. This leads to the equations 
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where A and Do are the gap function and diffusion con- 
stant respectively, rescaled by A , the zero temperature 
value of the gap. Finite temperature effects enter into 
the above equations through the spacing Alo = 2nT/ Ao 
of the rescaled Matsubara frequencies u> n . The Usadel 
equations are smooth in the frequency Lo n and hence the 
low temperature limit is expected to be reached when 
the temperature is about 10 times smaller than the gap 
A = 1.764fc B T c ~ 2.0 K in Al. We find by comparing 
our numerical results in the K-0 limit that an upper cut- 
off frequency of 8.0 K is sufficient for convergence. The 
equations are numerically solved self-consistently for A 



and g on a 2D grid. For 3D structures, we use a cylin- 
drically symmetric geometry to reduce it to an effective 
2D problem. The solutions on the 2D grid are calculated 
for different values of superconducting phase difference 
across the device and the current density is calculated 
using 



J(r) - -TrT^H^m^Hc^H) 



(2) 



where a is the conductivity of the metal film, e is the 
electron charge and the gradients are computed along 
the length of the bridge. 

The computed CPRs for an Al bridge with width 
W = 45 nm and three different lengths are shown in Fig. 
Q] Panels (a) and (b) compare bridges with 2D and 3D 
banks, respectively, at T = 0.15K. We take T c = 1.2K. 
The banks for each bridge are 600 nm in each lateral di- 
mension. For the 3D case, the banks also extend in the 
vertical direction as a hemi-cylinder of 600 nm diameter. 
For 3D structures, the CPR resembles a distorted sinu- 
soid for the shortest bridge length L = 15 nm~ 0.5£(0). 
As the bridge length exceeds £, the current varies more 
linearly with the phase and the value of the maximum 
current decreases. A similar trend is seen for 2D struc- 
tures of the same lateral dimension, but they exhibit a 
nearly linear CPR, even for L < £(0). To validate our nu- 
merical simulation, we compared our result for a 15 x 15 
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FIG. 2: Superconducting phase (color) as a function of posi- 
tion for an Al nanobridge with L = 75 nm and W — 45 nm. 
In (a) the banks are 2D with 600 nm lateral dimensions, and 
in (b) they are 3D as realized by a hemi-cylinder with 600 nm 
diameter. Contours of constant phase are also indicated. A 
total phase difference of 2 radians is imposed symmetrically 
across the bridge. In the 3D structure, most of the phase drop 
is across the constriction whereas in the 2D planar geometry, 
there is significant phase variation in the banks. 
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nm bridge contacted to ideal phase reservoirs with the 
analytical K-0 result and found quantitative agreement. 
When computing the current for 2D devices, the bridge 
cross-sectional area was chosen to match that of the cor- 
responding 3D devices by choosing the appropriate thick- 
ness for the 2D device. Note that the CPRs in general 
should be an odd function of S and also be 2n periodic 
which leads to the requirement that the current should 
be zero at S = ±mr [1J|. All the curves for the 2D de- 
vice and that for L — 240 nm for the 3D device have 
a non-zero value of the current at 5 = ir. However the 
full CPRs in these cases are multivalued and the other 
branch (suppressed in Fig. [1] for clarity) goes through 
zero at 5 = it. 

In order to understand the different shape of the CPR 
in the 2D and 3D geometries, we plot in Fig. [5] the su- 
perconducting phase as a function of position in a 75 nm 
long, 45 nm wide weak link junction — dimensions read- 
ily achieved with e-beam lithography. For these plots, 
a phase difference of 2 radians is imposed symmetrically 
at the ends of the banks (X = ±637.5 nm) and varia- 
tion with position is calculated. The phase is indicated 
by color and contour lines in Fig. O In the 3D case, 
the phase varies mostly in the bridge region and quickly 
heals to the imposed value in the banks within a few £(0). 
Thus, the banks act as good phase reservoirs and the 
structure has a nonlinear CPR. In contrast, for the 2D 
case, the phase evolves both in the bridge and the banks, 
with a logarithmic approach to the imposed value at the 
boundaries. Since the banks fail to act as phase reser- 
voirs, the entire structure resembles a superconducting 
wire with a linear CPR rather than a Josephson junc- 
tion. This is particularly a problem when multiple 2D 
bridges are incorporated in a more complicated circuit 
like a SQUID loop. There is significant overlap of the 
weak link phases resulting in reduced modulation depth 
[To| . Some of these limitations of 2D Dayem bridges have 
been discussed in [6j . Variable thickness bridges are thus 
a more appropriate replacement for tunnel junctions. 

We now consider both classical and quantum anhar- 
monic oscillators formed by capacitively shunting a 3D 
weak link junction. We neglect the small intrinsic ca- 
pacitance of the nanobridges Q and assume that the 
plasma frequency is determined by the additional shunt- 
ing capacitance only. For the classical case, we compute 
the oscillating junction voltage under a microwave drive. 
We consider a simple oscillator model where a weak link 
with generalized CPR 1(5) = I f(S) with max \f(S)\ = l 
is shunted in parallel with a capacitance Cs and a real 
impedance Rs which sets the quality factor of the os- 
cillator. The oscillator is driven with a time dependent 
current I(t) — Irf cos (u>dt) ■ The resulting equation of 
motion is 
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The equation of motion was solved numerically using a 
fourth order Runge-Kutta method for CPRs correspond- 



FIG. 3: Microwave frequency response of a nonlinear oscil- 
lator constructed using an Al nanobridge with L — 75 nm, 
W = 45 nm connected to (a) 3D and (b) 2D banks. The 
oscillator is driven with a single frequency u>d- The phase <f> 
of the steady state current oscillations, for a reflection setup, 
is plotted as a function of cud and normalized drive current 
Irf/Io- Increasing and decreasing current sweeps are inter- 
laced to highlight the hysteretic bistable region. Both devices 
exhibit pulling of the resonant frequency (yellow) to lower 
values with increasing drive current, typical of a driven non- 
linear oscillator with a softening potential. However, beyond 
the critical point P c (white dot), only the 3D device exhibits 
a well defined bistable region (striped) where the oscillations 
are confined to a single well of the periodic potential. 



ing to nanobridges with W = 45 nm, L — 75 nm con- 
tacted with 2D and 3D banks. The CPRs for these two 
devices differ in the maximum current Jo and the zero 
bias inductance Lj — <fo{dI(5)/d5)~ 1 \s=o. We chose 
values of Cs and Rs to obtain a small oscillation reso- 
nant frequency of 1.5 GHz and a quality factor Q = 50 for 
both samples — typical parameters in bifurcation ampli- 
fier circuits flo| . The results of the simulation are shown 
in Fig. [3^, b. The steady state oscillation phase, mea- 
sured in reflection, is plotted in color as a function of 
drive frequency u>d and normalized drive current Irf /Iq ■ 
Let us consider the variable thickness bridge first. At low 
drive amplitude Irf/Io < 0.01, one recovers the famil- 
iar linear resonance behavior as the phase evolves from 
180 to —180 degrees. The zero crossing of the phase 
(yellow) corresponds to the resonant frequency. As one 
drives the oscillator with larger drive current, the reso- 
nance frequency shifts to lower values as is expected for a 
softening potential. For currents beyond a critical value 
P c , indicated as a white dot in the figure, the oscillator 
becomes bistable and can coexist in two dynamical states 
(dashed region). ^This behavior is typical of driven nonlin- 
ear oscillators [l6| and has been experimentally verified 
in Josephson tunnel junction based oscillators [17| . This 
suggests that nanobridges with 3D banks can be used for 
potential applications in the field of bifurcation and para- 
metric amplification [18j. On the other hand, we do not 
observe a well defined bistable region in the microwave 
response of the Dayem bridge structure. This is because 
accessing the nonlinear regime requires larger phase ex- 
cursions and the phase S can hop between neighboring 
wells of the periodic potential [19j |. The critical point 
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FIG. 4: Quantum energy levels in a current biased, capaci- 
tively shunted L = 75 nm, W = 45 nm Al nanobridge with 
3D banks. For comparison, the potential and energy levels of 
a tunnel junction with the same critical current are shown as 
dashed lines. A current bias equal to 99% of the critical cur- 
rent is applied to obtain 3 quantum levels in the metastable 
potential well. The shunting capacitances are chosen to yield 
identical plasma frequencies. The level anharmonicity for the 
weak link junction (~ 3%) is only slightly smaller than that 
of the tunnel junction (~ 5%). 



can be pushed to lower drive currents but this requires 
significantly higher values of Q. 

Finally, we consider the possibility of constructing a 
quantum bit using weak link junctions. Since the critical 
current of these junctions tends to be large (~ 20^ A), 
it is natural to consider a phase qubit geometry in 
analogy with tunnel junction circuits with fiA critical 
currents [2(j. We will only present results for a 3D weak 
link junction and compare it with a Josephson tunnel 
junction since nanobridges with 2D banks have multival- 
ued CPRs which complicates calculations and the poten- 
tial operation of such a qubit. For the purpose of compar- 
ison, we also numerically solve the Schrodinger equation 
of capacitively shunted current biased Josephson junc- 
tions and find the quantum bound state energy levels. We 
choose a tunnel junction with a critical current equal to 
that of the weak link junction. An appropriate shunting 
capacitance is used to fix the zero bias plasma frequency 
of the tunnel junction device to ujpo/2ir = 20 GHz. A 
static current bias equal to 99% of Io is used to tilt the 
washboard potential to reduce the plasma frequenc y t o 
about 7.5 GHz, which is typical for phase qubits [2fj |. 
The 3D bridge with L = 75 nm and W = 45 nm is ca- 
pacitively shunted to achieve the same reduced plasma 



frequency under identical current bias. 

The potential wells and bound state energies are shown 
in Fig. |4j The tunnel junction data are shown as dashed 
lines. The overall shape of the potential for the two cases 
is quite similar. The positions of the minima occur at dif- 
ferent values of S, as indicated by two horizontal axes in 
the figure, since the maximum in the CPR of the tunnel 
junction is at 6 = n/2 and at S > tt/2 for the weak link 
junction. The lowest bound states have nearly the same 
energy with the subsequent higher levels differing by a 
few percent. For the tunnel junction, Wqi = 7.302 GHz 
and u>\2 — 6.948 GHz, corresponding to an anharmonic- 
ity of ~ 5%. For the weak link junction, wqi — 7.318 
GHz and u>\2 — 7.087 GHz, corresponding to an anhar- 
monicity of ~ 3%, which is only slightly smaller than 
that of the tunnel junction and suggests the plausibil- 
ity of weak link phase qubits. Moreover, it is possible 
to increase the anharmonicity by using a slightly higher 
plasma frequency. 

In conclusion, we have investigated the microwave 
transport properties and quantum energy level struc- 
ture of weak link Josephson junction oscillators, focus- 
ing in particular on optimizing anharmonicity. We com- 
puted the CPR for 2D Dayem and 3D variable thick- 
ness bridges by solving the Usadel equations for differ- 
ent bridge dimensions. The 3D structures exhibit strong 
nonlinearity — a consequence of the fact that the phase 
drop in these structures is confined mainly to the weak 
link. Under microwave drive, these junctions should ex- 
hibit bifurcation phenomenon for practical bridge dimen- 
sions and oscillator quality factor. Additionally, under 
static current bias, the quantum energy levels closely 
mimic ideal tunnel junction behavior, suggesting the pos- 
sibility of constructing phase qubits without any lossy 
oxide barriers, giving potentially enhanced coherence. 
Nanobridge junctions proposed here can also potentially 
be used to couple superconducting circuits to other solid 
state quantum systems such as molecular magnets and 
quantum dots. 
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